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1 Introduction 



The recent discovery of neutrino oscillations [IH, pointing towards non- vanishing neutrino 
masses, provides special motivation to investigate the question of fermion masses and mixing, 
at present one of the major riddles of particle physics. 

One of the striking features of the experimental evidence is the fact that large neutrino 
mixing is required at least in the 2-3 sector. This is to be contrasted to the situation in the 
quark sector, where the mixing is known to be small. In the search for a model which would 
naturally accommodate the experimental data, we find desirable to abide by the following 
principles: 

(i) We consider the standard model (SM) with the addition of three right-handed neu- 
trinos, but no extra Higgs doublets. The smallness of neutrino masses results then from the 
large right-handed neutrino masses through the seesaw mechanism [Q]. 

(ii) We will treat all the fundamental mass matrices on equal footing. In particular, we 
will assume that there is a weak-basis (WB) where, in leading order, the mass matrix of 
charged leptons and both the Dirac and right-handed Majorana neutrino mass matrices are 
all of the "democratic type", i.e., they are proportional to a matrix whose elements are, in 
leading approximation, all equal to unity. We will refer to this assumption as "extended 
democracy". 

(iii) We will assume that the breaking of extended democracy is small and it has the 
same pattern for all the fundamental mass matrices. 

In this letter we present a class of models which abide by the above principles and where 
the large mixing in the 2-3 sector results from the seesaw mechanism and requires a strong 
hierarchy in the right handed neutrino masses. 

It is worth emphasizing that assumption (ii) goes much beyond a simple choice of WB. 
Even if one assumes that both the neutrino Dirac and right-handed Majorana mass matrices 
have hierarchical eigenvalues, in general it does not follow that there is a WB where both 
the neutrino mass matrices as well as the charged lepton mass matrix are all, in leading 
approximation, proportional to the democratic matrix. If our framework is the correct one, 
it would mean that the appearance of large mixing in the leptonic sector originates in the 
seesaw mechanism. We recall that in the quark sector one can also obtain a good fit for the 
quark masses and mixings, assuming that both the up and down quark mass matrices are, to 
leading order, proportional to the democratic matrix, with a small perturbation generating 
the masses of the first two generations. The crucial new ingredient in the leptonic sector is 
the presence of the seesaw mechanism. 

Motivated by the recent discovery of neutrino oscillations, there have been in the lit- 
erature a large number of ansdtze for the structure of leptonic mass matrices . In most 
of them, a particular pattern is suggested directly for the effective left-handed neutrino 
mass matrix. A distinctive feature of the scheme we propose, is the fact that we suggest a 
universal pattern for all the fundamental leptonic mass matrices, which in our framework 
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are the charged lepton mass matrix, the Dirac neutrino mass matrix and the right-handed 
Majorana mass matrix. The structure of the effective left-handed neutrino mass matrix is 
then derived through the seesaw mechanism. 



2 General framework 



We consider the three generations SM, where three right-handed neutrino fields have been 
added, leading to the following charged lepton and neutrino mass terms: 

- Cmass = Ul iMi)ij IjR + PiL iMn)ij l^jR + ^^IrC {MR)ij Ujr + h.c. , (1) 

where the notation is obvious. Since the right-handed Majorana neutrino mass terms are 
SU{2)l X f/(l) invariant, Mr is naturally large, not being protected by the low energy 
symmetry. Following our general principles stated in the introduction, we will assume that 
all the matrices M^, Mr, and Mr are, to leading order, proportional to the democratic 
matrix A and, furthermore, that the breaking of the extended democracy has the same 
pattern for all the mass matrices. We thus write: 



Mfc = Cfe[A + Pfc], (2) 



where 
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with I ctfcl , I -C 1, so that all matrices are close to the democratic limit. As previously 
mentioned, the ansatz of Eqs. (0) and (0) is not just a choice of WB, together with the 
assumption of hierarchical masses for Mi, Md and Mr. Indeed if one assumes hierarchical 
masses, one can always choose, without loss of generality, a WB where, for example, both 
Ml and Mr are close to the democratic limit. However, once the vr basis is fixed, the Dirac 
neutrino mass matrix Mr, cannot in general be reduced to the quasi- democratic form by 
a choice of the basis. Thus, it is not possible in general to choose a WB where all the 
three matrices are in leading order proportional to A. Therefore, the nontrivial content of 
our ansatz of Eq. (|^) is the assumption that such a choice of WB is possible, implying an 
"alignment" of all three matrices in flavour space, and the suggestion that the breaking of 
the extended democracy has the same form for all three leptonic mass matrices. 

Following the hints of some Grand Unified Theories (GUTs), we consider the mass 
spectrum of M^) similar to the one of the up-type quarks. This will allow us to establish 
the relations between {a ,h ,c)d and the quark masses m^j, rric and mt- Taking into account 
that no Higgs triplets have been introduced, the effective mass matrix for the left-handed 
neutrinos is given by 

Meff = -Md Mr-' Md^ = -Ceff [A + Pd] Z [A + P^Y , (4) 
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where Cea = cd'^/cr and 



[ ciR + bR + ajibR -bji -an 

-bR bR I . (5) 

V - aR ttR 



It is convenient to define the dimensionless matrix Mq = —Mes/ces which can be written as 

Mo = AZA + AZPd + PdZA + PdZPd • (6) 

Due to the form of Z, the first term in Mq gives AZA = J2 Za A = A, the second and 

U J 

third terms vanish, while the fourth term reduces to 

/ \ 

PnZPo = X = Peff , (7) 
VO y) 

where x = a£)'^/aR and y = b£)'^/bR. The effective hght neutrino mass matrix can then be 
written as 

Meff = -Cefr[A + Peff]. (8) 

It is interesting to notice that this matrix has the same general form as the matrices Mi, 
Md and Mr, i.e. the seesaw mechanism preserves our ansatz. This is a remarkable feature 
of the scheme we propose in Eqs. (|^) and (^. 



3 Neutrino masses and mixing 

We shall choose the values of the parameters in the mass matrix Mcs so as to satisfy 
the experimental constraints on neutrino masses and mixings which can be summarized as 
follows. The Super-Kamiokande atmospheric neutrino data imply Am^g ~ (2 — 6) x 10"'^ 
eV^, sin^ 2623 > 0.84, and the combined data of the solar neutrino experiments lead to four 
domains of allowed values of Am^i and 612 corresponding to the four neutrino oscillation 
solutions to the solar neutrino problem - large mixing angle MSW (LMA), small mixing 
angle MSW (SMA), vacuum oscillations (VO) and low-Am^ (LOW) solutions @] []. For the 
remaining mixing angle 613, which determines the element of the lepton mixing matrix, 
only upper limits exist. The most stringent limit comes from the CHOOZ reactor neutrino 
experiment which together with the solar neutrino observations gives | sin 6*13! = l^esl < 
(0.22 - 0.14) for Aml2 = (2 - 6) x lO^^ eVl 

^Most recent Super-Kamiokande data disfavour the SMA and VO solutions at 95% c.l. However 
these solutions cannot yet be considered as ruled out and so we discuss them here along with the LMA and 
LOW solutions. 
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The hierarchical structure of the eigenvalues of the mass matrix of charged leptons Mi 
implies that it is very close to the democratic form, i.e. \ai\, \bi\ -C 1. The democratic mass 
matrix A can be diagonalized as F^AF = diag(0, 0, 3) with the real orthogonal matrix F: 
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The matrix Ui that diagonalizes Mi can therefore be written as 

U, = FW, 



(9) 



(10) 



where, due to the hierarchy \ai\ <^ ^ 1, the matrix W is close to the unit matrix. Next 
we will analyze two specific cases of the ansatz of Eq. (Q). 



3.1 The case of real mass matrices 

Let us consider that all the parameters and hk in Eqs. (0) and are real. It is instructive 
to consider first the limit when the matrix W in Eq. (|10D coincides with the unit matrix. 
The effective mass matrix of light neutrinos Meg, in the basis where the mass matrix of 
charged leptons has been diagonalized, is then obtained from Eq. (j^) through the rotation 
by the matrix F: M^s = F^M^gF. 

The first matrix on the r.h.s. of Eq. (|[) becomes diagonal upon this rotation. Therefore, 
if it dominates over the second matrix (i.e. if \y\ <C 1), all lepton mixing angles, including 
623, are small. This is phenomenologically unacceptable. Therefore we shall require that the 
second matrix in Eq. (H), i.e. P^s, either dominates or is of the same order as the first one. 
Since Pes is diagonal in the basis where Mi has an (almost) democratic form, in the basis 
where Mi has been diagonalized, Pes is non-diagonal and is diagonalized by the matrix F^. 
This means that, when Peff dominates in Eq. (^, the lepton mixing matrix U takes the form 
U ~ F'^, with F given by Eq. (|^). Therefore, in this case the mixing angle 612 responsible 
for the solar neutrino oscillations is Ou = 45°. One obtains also sin^ 2^23 = 8/9, which is 
within the range allowed by the Super-Kamiokande atmospheric neutrino data, and 613 = 0, 
in agreement with the CHOOZ limit. The value 6*12 = 45° is suitable for the VO and LOW 
solutions of the solar neutrino problem, but leads to slightly too high a value of sin^ 2^12 in 
the case the LMA solution which requires sin^ 26'i2 < 0.97 at 99% c.l. As we shall see, it is 
easy to satisfy this requirement if one takes into account a small contribution from the first 
matrix. A, in Eq. (P). For the SMA solution, the contribution from A in Eq. (||) should be 
comparable to that of Pes, otherwise 6*12 will be too large. 

We now proceed to analyze the effective mass matrix of the left-handed neutrinos. The 
parameters Ofc, bk and Ck in Eqs. (0) and are related to the masses of charged leptons, 
up - type quarks and heavy Majorana neutrinos through 

' rrir ' 2 mr ' 
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where we have taken into account the mass hierarchies present in the charged leptons and 
up-type quark sectors. The heavy right-handed neutrino masses are denoted by Mi, M2 
and M3. 

The effective mass matrix of in the basis where Mi has been diagonahzed takes the 
form 



Here 



eff = -CcSy 
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= -c^sy Mo . 
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3 Kmc/ Ml 



2 /mt\2 M2 

3 VmJ M3 



(13) 



The above discussed requirement, that the first matrix in Eq. (^ does not dominate over the 
second one, reduces to the condition \S\ < max{\e\, 1}. We also have to require \e\, \6\ -C 1 



in order to have the correct hierarchy Am 
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Amg ^ Am32 



The largest 



eigenvalue of the matrix Mq in Eq. (|12|) is then always close to unity. Thus, the value of 



Cesy (and so of M2) can be fixed by the requirement m| 



Am|2 



^^Itm, which gives 



m: 



4 X 10^° GeV 



(14) 



Using Eqs. (p!3|), it is then easy to find 



2 

Ml ~ - 



mt 



Am^ 



atm 



M. 



m: 



Am^ 



(15) 
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The masses of heavy Majorana neutrinos Mi, M2 and M3 are the only free parameters in our 
model and so all the neutrino masses and leptonic mixing angles can be expressed through 
them (or, equivalently, through qr, bR and M3 = 3|c/j|). 

Consider first the case \e\ <C \6\ -C 1 relevant for the SMA solution of the solar neutrino 
problem. In this case the eigenvalues of the neutrino mass matrix M^s in Eq. (|1^) are 



T fe 2^ £ ^ 5' 

jmi, m2, mg} ~ -c^sy <-, 3^ + 2' ^ + 3 



(16) 
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and the diagonalization of Mes results in the following lepton mixing matrix: 



U 
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(17) 



From Eqs. (|T6D and (0) we find 

£ ~ sin 26'i2 
g 

sin^ 2^23 - g 



2 

32 



^ Am 
2 
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2\ Ami2 



sin 6* 
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e3 



3v^ 



(19) 



Here we have used the fact that the smallness of |?7e3| implies sin^ 2^i2 — 4f/g^f/g2, sin^ 2^23 — 
4^^3^t3- a sample choice of the values of a^, and M3 that yields the SMA solution and 
the resulting neutrino parameters are shown in table 1. 

We shall now consider the case \b\ ;$ \e\ <^ 1 which, as we shall see, is relevant for 
the VO, LOW and LMA solutions of the solar neutrino problem. The eigenvalues of the 
neutrino mass matrix Meg in Eq. (|12D are 



{mi, m2, ms} ~ -c^sV 



S 



^ ^ 1 f 

3 9e' ^^3^9e' ^3 



(20) 



The diagonalization of M^s results in the following lepton mixing matrix: 



U 



( cos 6 
^sin^(l + |5) 



sin 6 

^COS0(1 + |5) 



3v^ 



V 



2sin^(l-|) 



cos 0(1 



(21) 



where 



tan 6' 



2 5 252 

3£ ^ 9^2 



:22) 



In the limit \5\ <C \e\, tan^ 1 and the mixing matrix U of Eq. (|21|) becomes the 
matrix (up to the trivial sign changes due to the rephasing of the neutrino fields), as it 
should. Notice that for 5 7^ the element f/e3 of the lepton mixing matrix is no longer zero. 

From Eqs. (|20|) - ( p2|) it is easy to find that, to leading order in [/e3, 

tan 6112 = tan 61 , (23) 
V2 tan 6^12 — 



tan6'i2x— §^ , 



5 



tan0i2x— §^ , 



(24) 



^We use the parametrization of the leptonic mixing matrix U which coincides with the standard 
paramctrization of the quark mixing matrix . 
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which substitute for Eq. (pTSD, whereas Eq. ([T9|) remains vahd in this case. 

Eqs. ( p^ - ( p^ allow one to find the values of a^j, bn and M3 necessary to obtain the 
relevant solutions of the solar neutrino problem. Sample choices of these parameters, along 
with the resulting values of neutrino masses and lepton mixing parameters, are shown in 
table 1. For example, if one chooses = 1.9 x 10"^°, 6r = 8.0 x 10"^ and M3 = 2 x 10^*^ 
GeV, one obtains e = 5 x 10^^, 5 = 2.3 x 10^^ which leads to the LOW solution of 
the solar neutrino problem. In particular, the solar neutrino mixing is nearly maximal, 
2 2^19 = 0.999. 



sm 



The value 612 — 45° is also appropriate for the VO solution, but for this solution the 
best fit corresponds to sin^ 2^12 — 0.7. Therefore, we require 5 ~ £ to have a smaller value 
of sin^ 2^12 than in the case of the LOW solution. Notice that, for the VO solution, we also 
have to require 5 ~ e in order to avoid unacceptably large values of M3 (see discussion in 
sec. 4). A choice for the parameters qr, Br and M3, very close to the best fit VO solution 
for Am^i and sin^2^i2, is shown in the last column of table 1. 

The LMA solution also requires 5 ~ e in order for 612 not to be too close to 45°. A 
possible choice of the values of qr, Br and M3 and the resulting neutrino parameters are 
shown in table 1. Notice that in this case sin^ 2^23 — 0.91 which is slightly larger than the 
value 8/9 obtained for the LOW and VO solutions. These values are still slightly smaller 
than the ones obtained for the SMA solution. From Eq. (|19D, it is clear that this is due to 
the fact that, from all the solutions, the SMA requires the largest value of 6. 

We shall now take into account that the mass matrix of charged leptons is not exactly 
of the democratic form. Thus, the matrix W in Eq. ( |lOl) deviates slightly from the unit 
matrix due to nonzero values of rrie and . Diagonalization of Mi = F'^MiF gives 



W 



( 



V 



^/3n 



3 r^e 
2 rrir 
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2 me \ 

3 nir 



V2v 



(25) 



The lepton mixing matrix is now obtained from Eq. (|2l| 



Q replacing U W^U. 

This modifies the lepton mixing parameters. In the case \S\ < \e\ <^ 1 relevant for the VO, 
LOW and LMA solutions of the solar neutrino problem one finds 



or Eq. 

< 



sin^ 2^12 = sin^ 26 
sin^ 2^23 = l(l + ^6 
sin 6^13 



4 m 
1 - - 



cos 20 



9 

Ue3 



2 
3 

e6 



3 1 — cos 29 
1 + !^(1_35) 



x/2 



nif. 



3^2 3 



In the limit me, 



the corresponding expressions of Eqs. 



(26) 
(27) 
(28) 

and (E3D are recovered. 
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Analogously, in the case \e\ ^ \S\ ^ 1 relevant for the SMA solution of the solar neutrino 
problem, one obtains from Eq. (|r^) and Eq. 



sm' 



'26 



12 



3e 
46 



1 me 
3 m,, 



(29) 



whereas sin^ 26*23 and Ues are again given by Eqs. (^) and (^Sj) . 

The numerical values of the lepton mixing parameters with contributions from the 
charged lepton masses included are given in table 2. It is interesting to notice that these 
contributions tend to increase sin^ 2^23, bringing it closer to the Super-Kamiokande best fit 
value. They also increase significantly the value of the mixing parameter \Ue3\- Here, these 
contributions are dominant for all the solutions of the solar neutrino problem except, for 
LMA, in which case they constitute about 2/3 of |t/e3|. 



3.2 Universal strength of Yukawa couplings 

We shall now consider a special case of complex parameters and 6^ , which are of the form 

ak = e"'" - 1 , hk= e'"" - 1 , (30) 

with real ak and (3k- In this case all matrix elements of the matrices in Eq. @ have 
the same absolute values, i.e. this is the special case of the so called universal strength of 
Yukawa couplings (USY) p. 

The effective mass matrix of ul can again be written in the form of Eq. (^ with Pefr 
given by Eq. (^), where x and y are now complex parameters: 

al 2sm'{ao/2) 2sin^(/3^/2) 

/ / \ 5 y 7 //0/\ \ / 

ttR sin(ai?/2) bR sm{PR/2) 

with p = f + - ^ and 7 = f + - ^. 

The phases ak, (3k and the parameters Ck satisfy, to leading order, the relations of 
Eq. (pi]) with the substitution ak ^ ak, bk ^ (3k- 

Next, we define the complex parameters e' = x/y and 6' = 3/y, and denote e = \e'\, 
S = Since Mi is almost democratic, it can be approximately diagonalized by the matrix 
F. In this approximation,^ the effective left-handed neutrino mass matrix Meg, in the 
basis where Mi has been diagonalized, is given by Eq. (|r^) with e and 6 substituted by 
e' and 6'. The leptonic mixing matrix then coincides with the unitary matrix U that 
diagonalizes Mes- To find this matrix it is convenient to define the Hermitian matrix 
Hes = Meff^iff = \ces yl"^ MqMq , which can be diagonalized through V^^HqV. The matrices 
U and V are then related by f/ = V*K, where if is a diagonal matrix of phases. 

^This approximation amounts to neglecting the small terms of the order of nie/m^, me/rrir and m^/mr 
in the unitary matrix W defined in (10), setting it to the unit matrix. 
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Further simplification can be achieved by noticing that the phases p and 7 are very close 
to 7r/2, which gives e' ~ e and 5' ~ —i5. The matrix Mq is then obtained from Eq. ( [T^ ) 
substituting 5 — > —i5. All its matrix elements except (Mo)33 are real. 

We shall first consider the case \e\ <C \5\ -C 1 which is relevant for the SMA solution of 
the solar neutrino problem. In this case, the eigenvalues of ifes are 



{ml, ml, ml) ~ |ccfri/p — , -5^ + -e^, 1 + 



9 



4 



9 



(32) 



where the mj's are the effective left-handed neutrino masses. The diagonalization of Hq 
leads to the following lepton mixing matrix U: 



U 



'AS 



- (I 

4 <5 '^-^ 



l5) 



V3 



''3V2 ^ 
il-^6) 



\ 



3 e 
8 <5 



1 
V3 



(33) 



As in the case of real mass matrices, in the limit \e\ <^ \S\ <^ 1 the parameters e and S are 
related to physical observables through Eq. (13). For the mixing angles one obtains 

— gin ^V- ^ I - 

9 ' ' 



sin^^ 



12 



. 2n 



131 



(34) 



Choosing for the SMA solution the same input parameters as in the case of real mass 
matrices, we get the results similar to those obtained in the previous subsection (see table 3). 

Next, consider the case \6\ ^ \e\ <^ 1, relevant for the LMA, LOW and VO solutions. 
We get the following results for the eigenvalues of Hq: 

.f. 5^ 5^ ' 



3^2 



1 + 
9 



(35) 



Introducing k, = 6/3e one finds for the mixing matrix U: 



/ ^(1 -2m-4/t4) ^ (1 - 2m - 4/t2 + + 12k^) \ 



U 



V2 



2k'^ + 2k^) 

Notice that k < 1/3. The mixing angles are 
sin 2^12 ~ 1 - 16k^ + 32/t6 , 



-^(1^2k^-Qk') 



-^{1 + 2k^-6k^) 



3V2 
1 

V3 



(36) 



sin^ 2^23 



9 



sin 6 



13 



e6 

3V2' 



(37) 



In table 3 we give the results of the numerical diagonalization of the effective neutrino mass 
matrix (with no approximations made) and compare them with those obtained using the 
approximate analytic expressions. 



9 



The results for the USY case, presented so far, were obtained neglecting the corrections 
of the order mg/m^, me/mr and m^/rrir coming from the diagonalization of Mi, i.e. by 
replacing W by the unit matrix. Consider now these small corrections. The lepton mixing 
matrix will be U' = W'^U . The unitary matrix W is (up to a diagonal phase transformation) 



the matrix that diagonalizes H'l 
obtained the following result: 



F^MiM}F as W^H[W = diag(m; 



65 ^fl^ 



m: 



'). We have 



W 



( 



\ 



1 + 

2 m-r 



—I 



2 me ^ 
3 



4 m-r 4 JTtr 



3 nie 
2 m-r 



' \/2mT 



(St 



J 



It is easy to show that, unlike in the case of real mass matrices, the corrections coming 
from W ^ 1 give negligible contributions to the mixing angles 612 and 623- This can also 
be seen by comparing tables 3 and 4, which give the results of the numerical simulations 
performed with and without contributions of charged lepton masses. These corrections, 
however, constitute the main contributions to |sin6'i3|. Indeed, taking them into account 
one finds 



sin^ 



13 




6e nip 6"^ e"^ 



+ 7r— + 



(39) 



which is valid for all the solutions (SMA, LMA, LOW and VO). In the limit mg/m^, me/m,- — »• 
one recovers the result given in Eqs. (^4|) and (^7|). However, for the realistic values of 



parameters, the contributions of the order of me/rrir and {me/m^Y ^"^^ always important. 

We shall now discuss briefly the leptonic CP violation. In general, in the case of three 
light Majorana neutrinos there is a Dirac-type CP violating phase 5cp and two additional 
Majorana-type phases. The latter cannot be observed in neutrino oscillation experiments 
and we shall not discuss them. The phase (5cp can be found from the invariant CP violating 
parameter \J'\ 



llm 



and the values of the mixing angles. The values of 
J and (5cp, calculated numerically for all four solutions of the solar neutrino problem, are 
presented in tables 3 and 4 0. 



3.3 Renormalization group effects 

So far, we have analysed the phenomenological implications of an ansatz for the charged 
lepton and neutrino mass matrices, based on the hypothesis of extended democracy. In 
this section, we study the behaviour of the ansatz under the renormalization group. This 
analysis is especially important due to the fact that if the observed pattern of fermion masses 

'^The values of J obtained from the approximate expressions (^3|) and (^6|) are zero. In order to get 
non-zero results as presented in table 3 one has to consider higher order terms in the mixing matrices which 
do not contribute significantly to the mixing angles. 
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and mixings reflects some flavour symmetry of the lagrangian, it is natural to assume that 
this symmetry is manifest at a high energy scale and analyse its predictions at low energies 
by studying its evolution under the renormalization group. In the SM context, the lowest 
dimension operator that can generate a Majorana mass term for the left-handed neutrinos 
is uniquely given by ^ 

--Ku^uHH + h.c (40) 



where k is proportional to Mq in Eq. (|T2D, and H is the neutral component of the usual SM 
Higgs doublet. Assuming that n is defined at Mr, we will study the stability of the implied 
pattern of neutrino masses and mixings by analysing the renormalization group equation 



(RGE) for the operator k which can be written, at one loop level, as |10]: 



lOTT — 

dt 



2A + 6Y^^ - 3g. 



K 



K 



Yt Ye + (YtYe)^ 



K 



(41) 



where t = log A, and g2, X,Yt,Ye, A are the SU{2) gauge coupling, the quartic Higgs 
coupling, the top Yukawa coupling, the charged leptons Yukawa couplings matrix and the 
scale at which k is evaluated, respectively. This equation allows us to relate the effective 
mass matrix at Mr with the one at low energies, say at the scale of the Z boson mass Mz, 



in the following way ||11|| : 



Mo{Mz) ~ A diag(l, 1,1 + v) Mo{Mr) diag(l, 1, 1 + r^) 



(42) 



where Mq{Mz), Mq{Mr) are the matrices Mq at the scales of Mz and Mr, respectively, and 
rj is approximately given by: 



V 



Yi 



327r2 



log 



A(Mj 
Mz 



1 /rrij 
327r2 \~ 



\ HMn) 
log 



M. 



(43) 



where rrir is the r mass and v denotes the vacuum expectation value of the neutral Higgs 
doublet component It is straightforward to see that the transformation defined in Eq. ( |1^ ) 
corresponds to multiplying the third row and the third column of Mo(Mr) by the parameter 



1 + 77. The overall factor A in Eq. (|42|) comes from the terms in the first square bracket 
on the right-hand side of Eq. (^) and it will affect the values of the neutrino masses but 
not the structure of the mass matrix. Taking into account the form of Mo(Mr;) given in 
Eq. (0), we get, neglecting second order terms: 



Mo(Mz) (X Mo{Mr) + 



f 













1 



2r] 



J 



(44) 



In the previous sections, we saw that our results implied a large hierarchy between the 
right-handed neutrino masses which can lie between 10® GeV and 10^^ GeV, corresponding 



^ We consider 



1.777 GeV, Mz = 91.187 GeV and v = 174 GeV. 
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to the lightest and heaviest j^Ft, respectively (see tables 1-4). Using this range for the scale 
A{Mr), we get 3.1 x 10"^ < < 1.2 x 10"^ Since 10"^ < e,S < 10"^ (see tables 1 - 4), it 
can be readily verified that even for rj = 1.2 x 10"^ the effect of the RGE on structure of 
the effective neutrino mass matrix is negligible. The point is that although the overall scale 
factor of Meff may run significantly, the structure of this matrix is quite stable. Therefore, 
the results we have obtained at low energies will still be valid if we impose the hypothesis 
of extended democracy, with universal breaking for all the leptonic mass matrices, at a high 
energy scale. 



4 Discussion 



We have proposed a model for lepton masses and mixing with a high predictive power. With 
just three free parameters (the masses of heavy right-handed neutrinos Mi, M2 and M3 or, 
equivalently, a^, and M3) we predict seven physical quantities - the masses of light 
neutrinos mi, m2, m^, the mixing angles 6*12, 6^23 and 6'i3, and the CP- violating phase 5cv- 
Depending on the values of the input parameters, all four solutions of the solar neutrino 
problem (SMA, LMA, LOW and VO) can be obtained. Representative choices for the input 
parameters and the resulting values of light neutrino masses, lepton mixing angles and CP- 
violating parameters are given in tables 1-4. We have performed both exact numerical 
and approximate analytic diagonalizations of the light neutrino effective mass matrix. As 
follows from tables 1-4, our analytic expressions give a very accurate approximation to the 
exact results. We have also demonstrated that the structure of the effective mass matrix of 
light neutrinos is stable against the RGE evolution. 

A characteristic feature of our scenario is a large hierarchy of all neutrino masses, in- 
cluding the masses of heavy Majorana neutrinos [jl2|. This is because we assume all the 



fundamental lepton mass matrices to have a nearly democratic form. It is interesting to no- 
tice that, for all the solutions of the solar neutrino problem, the values of M2, i.e. the mass 
of the second heavy Majorana neutrino, are nearly the same. This stems from the fact that 
M2, approximately given by Eq. ([T^), is related to Am^j^, and is practically independent 
from Ami, and 6'i2. 

Physically, the mass M3 of the heaviest right handed Majorana neutrino must not exceed 
Mpi = 1.2 X 10^^ GeV. Eq. (||) then gives a lower bound on \5\: \5\ > 4.2 x 10"^ However, 
for real lepton mass matrices, in the case of the LMA, LOW and VO solutions, it follows 
from the second equation in (|2^) that \5\ decreases with tan ^12 — >■ 1. Therefore values of 
^12 too close to 45° are not allowed in our scenario. For the LMA and LOW solutions, 
this restriction is not severe (the values of sin^ 2^12 as large as 1 — 4 ■ 10~^ for LMA and 
1 — 3 ■ 10^^ for LOW are allowed), but for the VO solution sin^ 26^12 must not exceed 0.967. 
Interestingly, in this case, the value of sin^ 26'i2 giving the best fit of the solar neutrino data, 
is not close to 1. As follows from the first equality in (^), in the complex USY case sin^ 26'i2 
can be very close to unity. 
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It should be noticed that for LOW and VO solutions M3 ~ 10^^ GeV (see tables 1-4), 
i.e. is close to the reduced Planck scale which is presumably the string scale. In these cases 
our results might be affected by new physics at this scale. 

It can be seen, if one compares tables 1 and 2 (and Eqs. ( ^^ ) and ([27|)), that in the case 
of real neutrino mass matrices the corrections of order m^/rriT-, coming from the deviation 
of the charged lepton mass matrix from the exact democratic form, increase the values of 
sin^ 2^23, bringing them closer to the Super-Kamiokande best-fit value sin^ 26*23 = 1. At the 
same time, in the case of complex parameters the corrections to sin^ 2^23 are of the order 
(m^/m,-)^, i.e. they are negligible (compare tables 3 and 4). This comes about because the 
terms of the order m^/mT- in the matrix W in Eq. (|38| ) are purely imaginary, unlike those 
inEq. (g). 

The corrections to sin^ 26*12 due to nonzero rrif. and are small (of the order me/m^) 
in the case of real lepton mass matrices and totally negligible in the case of complex USY 
matrices. In contrast to this, the corresponding contribution to jf/esl = Isin^isl, though 
of the order mg/m^, is dominant (compare tables 1 and 2, 3 and 4 and also Eqs. ([19| ) 
and (p8|)). For the SMA, LOW and VO solutions of the solar neutrino problem we find 
|?7e3| ~ (\/2/3)(me/m^) ~ 2.3 x 10^^, whereas for the LMA solution it is slightly larger. 
Unfortunately, these values are too small to be experimentally probed in currently planned 
long-baseline neutrino oscillation experiments. 

The values for \Uez\ that we have found are different from the predictions of Ref. []T3| 
obtained under the assumption of no fine tuning between the elements (Meg) 12 and (Meff)i3 
of the neutrino mass matrix in the basis where the mass matrix of charged leptons has 
been diagonalized. The reason for this is that, in our case, there is an approximate equality 
(Meff)i2 sin ^23 + (-^eff)i3 COS 6*23 — 0, which is exactly the kind of relation which was excluded 
from the consideration in [jl3[ . This relation can be traced back to an approximate symmetry 



underlying Cmass in Eqs. (|l]) - (||). Thus, our scheme provides an example of the case in 
which the predictions of [jl3| do not apply. 

In the case of complex lepton mass matrices, we predict relatively large values for the 
CP-violating phase ^cp in the case of SMA and VO solutions, and small values in the 
case of LMA and LOW solutions. The contributions due to nonzero me and are very 
important in this case - they increase the CP- violating parameter by 2 - 6 orders of 
magnitude for the SMA, LOW and VO solutions and decrease it by 2 orders of magnitude 
for the LMA solution. Unfortunately, CP-violating effects in neutrino oscillations cannot be 
experimentally probed in our scheme because of the smallness of which is mainly due 
to the smallness of the mixing angle ^13. 

In conclusion, we have suggested a simple structure for the leptonic mass matrices, 
with the remarkable feature that a simple explanation is provided for the large mixing in 
the leptonic sector, in contrast with the quark sector. It is well known that, in the quark 
sector, one may obtain the correct mass spectrum and mixing pattern starting with 
a democratic matrix for the up and down quarks and adding a small perturbation which 



13 



generates the masses of the two hght generations, as well as the small mixing present in 
the Cabibbo-Kobayashi-Maskawa matrix. In the scheme we have proposed, all leptonic 
mass matrices are treated in an entirely analogous way, i.e. they are, in leading order, all 
proportional to the democratic matrix, with a small universal breaking of democracy. The 
large mixing in the leptonic sector results from the seesaw mechanism, which is the crucial 
new ingredient, only present in the leptonic sector. 

We are grateful to F. Feruglio and M.N. Rebelo for useful discussions and to the referees 
of the paper for constructive suggestions. The work of E.A. and F.R.J, has been supported 
by Fundagao para a Ciencia e a Tecnologia under the grants PRAXIS XXI/BCC/16414/98 
and PRAXIS XXI/BD/18219/98, respectively. 
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LMA 


SMA 


LOW 


vo 


H 


aR 


1.5x10"' 


8.5x10"* 


1.9x10"'" 


4.1x10"' 


INPU 


bR 


1.3x10"^ 


2.2x10"^ 


8.0x10"* 


3.9x10"* 


M3 (GeV) 


1.3x10'* 


7.6x10'^ 


2.0x10'* 


4.3x10'* 




E 


1.04x10"' 


3.11x10"' 


5.06x10"' 


1.14x10"^ 




5 


3.72x10"^ 


6.30x10"^ 


2.29x10"^ 


1.12x10"^ 






Numerical 


Appr. 


Numerical 


Appr. 


Numerical 


Appr. 


Numerical 


Appr. 




Ml (GeV) 


3.2x10* 


3.2x10* 


1.1x10* 


1.0x10* 


6.3x10' 


6.3x10' 


2.9x10' 


2.8x10' 




M2 (GeV) 


3.8x10'" 


3.7x10'" 


3.7x10'" 


3.6x10'" 


3.6x10'" 


3.6x10'" 


3.8x10'" 


3.7x10'" 


H 


Am^i (eV^) 


5.36x10' 


5.38x10-' 


7.25x10-' 


7.54x10-* 


1.15x10-' 


1.15x10' 


1.02x10'" 


1.03x10'" 


OUTPl 


Ain^32 (eV^) 


3.94x10"^ 


3.94x10-' 


4.15x10-' 


4.14x10' 


4.35x10' 


4.35x10-' 


3.96x10' 


3.96x10' 


Sin ZO12 


0.95 


0.95 


5.69x10' 


5.48x10' 


0.999 


0.999 


0.70 


0.73 




sin^ 2023 


0.91 


0.91 


0.92 


0.93 


0.889 


0.889 


0.889 


0.889 




mi(eV) 


6.76x10" 


6.82x10-" 


9.44x10-' 


9.81x10' 


4.96x10-* 


4.96x10* 


1.64x10* 


1.58x10* 




m2(eV) 


7.35x10-^ 


7.37x10-' 


2.69x10-' 


2.75x10-' 


3.39x10-" 


3.39x10-" 


1.02x10-' 


1.03x10-' 




mj (eV) 


6.32x10-^ 


6.32x10-^ 


6.45x10-^ 


6.44x10-^ 


6.59x10-^ 


6.59x10-^ 


6.29x10-^ 


6.29x10-^ 




lUeal 


1.02x10"' 


9.14x10-" 


4.63x10-' 


4.62x10-' 


2.75x10-' 


2.73x10-' 


3.01x10-' 


3.01x10' 



Table 1 - Results for exact numerical and approximate diagonalizations of M^^i in the case of real lepton mass 
matrices, with sample choices of input parameters leading to the different solutions of the solar neutrino 
problem, not including the corrections due to nonzero and . 







LMA 


SMA 


LOW 


VO 


H 


aR 


1.5x10"' 


8.5x10"* 


1.9x10"'" 


4.1x10"' 


INPU 


bR 


1.3x10"^ 


2.2x10"^ 


8.0x10"* 


3.9x10"* 


M3(GeV) 


1.3x10'* 


7.6x10'^ 


2.0x10'* 


4.3x10'* 




E 


1.04x10"' 


3.11x10"' 


5.06x10"' 


1.14x10"^ 


3UTPUT 


s 


3.72.x ICr 


6.30.\10"" 


2.29x10""' 


1.12xKr^ 




Numerical 


Appr. 


Numerical 


Appr. 


Numerical 


Appr. 


Numerical 


Appr. 


sin^ 2012 


0.95 


0.94 


5.1x10' 


5.0x10' 


0.999 


0.999 


0.70 


0.72 




sin^2e23 


0.95 


0.96 


0.96 


0.97 


0.94 


0.94 


0.94 


0.94 




lUeal 


3.27x10-' 


3.17x10-' 


2.27x10-' 


2.28x10-' 


2.29x10' 


2.28x10-' 


2.29x10' 


2.28x10' 



Table 2 - Changes to table 1 if one includes the corrections due to nonzero and . 
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LMA 


SMA 


LOW 


vo 


INPUT 


Or 


3.8x10"' 


8.5x10"* 


1.9x10"'" 


4.1x10"' 


Pr 


2.5x10"' 


2.2x10"' 


8.0x10"* 


3.9x10"* 


M3(GeV) 


6.5x10" 


7.6x10" 


2.0x10'* 


4.3x10'* 


OUTPUT 


E 


7.91x10"^ 


3.11x10"' 


5.06x10"' 


1.14x10"^ 


5 


7.16x10"' 


6.30x10"' 


2.29x10"^ 


1.12x10"^ 


M, (GcV) 


4.1x10" 


l.lxio' 


6.3x10" 


2.9x10' 


M2(GeV) 


3.6x10'° 


3.7x10'° 


3.6x10'° 


3.8x10'° 




Numerical 


Appr. 


Numerical 


Appr. 


Numerical 


Appr. 


Numerical 


Appr. 


Am^i (eV^) 


3.26x10"' 


3.24x10-' 


7.05x10-*^ 


7.04x10-* 


1.11x10-' 


1.11x10' 


6.55x10" 


6.61x10" 


Aiii^2 (eV^) 


4.18x10-^ 


4.22x10-^ 


3.97x10-' 


3.99x10-' 


4.34x10' 


4.34x10' 


3.95x10' 


3.95x10' 


sin^ 2012 


0.91 


0.89 


5.48x10-' 


5.48x10-' 


1-8.3x10-' 


1-8.3x10' 


0.89 


0.86 


sin^ 2023 


0.89 


0.89 


0.89 


0.89 


0.89 


0.89 


0.89 


0.89 


nil (eV) 


1.35x10-^ 


1.32x10-^ 


9.78x10-' 


9.81x10-' 


5.03x10-"^ 


5.03x10-"* 


2.02x10-" 


1.93x10-" 


(eV) 


5.87x10-^ 


5.85x10-^ 


2.66x10-' 


2.66x10-' 


3.34x10-" 


3.34x10-" 


8.34x10-" 


8.36x10-" 


mj (eV) 


6.49x10-^ 


6.50x10-^ 


6.31x10-^ 


6.31x10-^ 


6.59x10-^ 


6.59x10-^ 


6.29x10-^ 


6.29x10-^ 


|Ue3| 


1.24x10-^ 


1.33x10-^ 


4.61x10-' 


4.62x10-' 


2.72x10-' 


2.73x10-' 


3.01x10-' 


3.01x10-' 


J 


-2.35x10-'' 


-5.92x10-* 


-6.41x10-* 


-4.87x10-'° 


IScpIO 


57.5 


4.2 


88.3 


46.7 



Table 3 - Results for exact numerical and approximate diagonalizations of M^ff in tlie case of complex USY 
lepton mass matrices, with sample choices of input parameters leading to the different solutions of the solar 
neutrino problem, not including the corrections due to nonzero and . The CP-violating parameters J" and 
I 5cp I , were also computed. 







LMA 


SMA 


LOW 


VO 


H 




3.8x10"' 


8.5x10"* 


1.9x10"'° 


4.1x10"' 


INPU 


Pr 


2.5x10"' 


2.2x10"' 


8.0x10"* 


3.9x10"* 


MaCGeV) 


6.5x10" 


7.6x10" 


2.0x10'* 


4.3x10'* 




sin' 2012 


0.91 


5.49x10' 


1-2.8x10' 


0.89 




sin' 2O23 


0.89 


0. 89 


0.89 


0.89 


OUTPU 


J 


3.44x10-" 


3.99x10' 


-7.75x10-" 


2.59x10-* 


|8cp|0 


0.3 


81.0 


0.8 


30.6 


|Ue3| 


Numerical 


2.83x10-' 


2.31x10-' 


2.29x10-' 


2.29x10-' 




Approximate 


2.69x10' 


2.28x10-' 


2.28x10-' 


2.28x10' 



Table 4 - Changes to table 3 if one includes the corrections due to nonzero and . 
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